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— A  basic  inverse  problem  in  underwater  acoustics  is  the  determination  of  the  structure  of 
the  sea  bed  from  a  limited  knowledge  of  its  reflection  coefficient.  For  many  applications,  an 
adequate  model  for  studying  the  acoustic  interaction  is  provided  by  the  scattering  of  plane 
waves  by  a  layered  liquid  medium.  In  contrast  to  formally  exact  solutions  to  this  inverse 
scattering  problem,  Candel,  et  al.,  [  J.  Sound  Vib.  68,  571-595  (1980)]  developed  an  ap¬ 
proximate  scheme  which  is  readily  implemented  numerically.  Under  the  assumption  of 
nonabsorptive  media,  the  reconstruction  of  both  sediment  density  and  sound-speed  profiles 
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ys  effected  by  the  niunerical  integration  of  a  system  of  four  first-order  differential  equations 
requiring  impulse  response  data  for  two  distinct  grazing  angles.  In  this  paper,  the  effect  of 
additive  noise  on  the  inversion  algorithm  is  examined  for  synthetic  reflection  data  generated 
for  a  geoacoustic  model  of  deep-sea  sediments  representative  of  the  Hatteras  Abyssal  Plain. 
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THE  EFFECT  OF  NOISY  REFLECTION  DATA  ON  AN 
INVERSE  METHOD  FOR  DETERMINING  THE  STRUCTURE  OF  A 
LAYERED  OCEAN  BOTTOM 


At  low  frequencies,  sound  waves  propagating  within  the  sea 
can  readily  penetrate  the  sea  bed  to  considerable  depths. 
Hence,  proper  interpretation  of  acoustic  signals  which  pene¬ 
trate  the  sea  bottom  requires  a  knowledge  of  the  density  and 
sound-speed  profiles  below  the  sea  floor.  A  basic  problem  is 
the  determination  of  the  structure  of  the  sea  bed  from  a 
limited  knowledge  of  the  reflected  waves.  For  many  appli¬ 
cations.  an  adequate  model  for  studying  this  inverse  problem  is 
the  scattering  of  plane  acoustic  waves  by  a  one-dimensional 
inhomogeneous  medium. 


THE  EFFECT  OF  NOISY  REFLECTION  DATA  ON  AN 
INVERSE  METHOD  FOR  DETERMINING  THE 
STRUCTURE  OF  A  LAYERED  OCEAN  BOTTOM 

BACKGROUND 

•  BURRIDGE  (1980) 

•  NEWTON  (1981) 

•  CAN  DEL,  DEFILLIPI  &  LAUNAY  (1980) 

•  THOMSON  (1983) 

Viewgraph  1 

Formal  methods  for  solving  one-dimensional  inverse  scat¬ 
tering  problems  are  reviewed  by  Burridge  (1980)  and  Newton 
(1981).  Candel.  Defillipi  and  Launay  (1980)  proposed  an 
approximate  scheme  which  is  readily  implemented  numerically. 
An  important  feature  of  the  Candel.  et  al.,  approach  is  the 
separate  recovery  of  both  density  and  sound-speed  profiles  from 
two  band-limited  impulse  responses  associated  with  different 
grazing  angles.  The  inversion  of  noise-free  reflection  data 
using  their  algorithm  has  been  presented  elsewhere  by  Thomson 
(1983).  However,  measured  sediment  reflection  data  will  always 
contain  some  level  of  noise  because  of  ambient  noise  in  the 
ocean  and  data  processing  limitations.  Today  we  will  present 
the  results  of  a  study  of  the  effects  of  inverting  simulated 
reflection  data  containing  various  levels  of  Gaussian  noise. 
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Model  for  One-dimensional  Inverse  Method 


z-0 


Z=H 


•  PLANE  WAVES 

•  HORIZONTAL  LAYERING 

•  LIQUID  LAYERS 

•  ZERO  ABSORPTION 


Viewgraph  2 

In  the  mathematical  model  to  be  considered,  the  hori¬ 
zontally  layered  region  0<z<H  is  occupied  by  an  inhomogeneous 
liquid  having  arbitrary  variation  of  density  p(z)  and 
sound-speed  c(z).  The  half-spaces  z  <  0  and  z  >  H  are 
homogeneous.  All  regions  are  assumed  to  be  non-absorbing. 

A  plane  acoustic  wave,  initially  propagating  downward  at 
the  grazing  angle  0O  encounters  the  inhomogeneous  medium 
at  z  *  0.  Within  0  <  z  <  H.  the  incident  wave  is  partially 
reflected  and  a  reflected  wave  is  returned  to  the  region  z<0. 
In  z  >H.  only  a  transmitted  wave  continues  to  propagate. 
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Frequency  Response  of  the  Reflection  Coefficient  as 
a  Non-linear  Fourier  Transform 

H  z 

•  R(f,8)  «  -  /  g(z)  exp[2i/kz  (s)ds)  dz 


1_  dYz 
2YZ  dz 


•  g  =  — — z  = 

®  At/ 


2  dz 


•  Y2  as(kz  Ik)  I  (pc) 

2  2  V* 

•  kz  =  (k-k*  )  =  k  sin  8 


Viewgraph  3 

As  a  result  of  applying  the  forvard  scattering  approxi¬ 
mation  to  a  local  wave  representation  of  the  field.  Candel.  et 
al..  derived  an  approximate  form  for  the  reflection  coefficient 
R  as  a  function  of  the  frequency  f  and  grazing  angle  Q.  In 
the  upper  line,  observe  that  R  appears  as  a  non-linear  Fourier 
transform.  The  kernel  of  the  transform,  g.  is  one-half  the 
logarithmic  derivative  of  "longitudinal  admittance",  Yz.  Yz 
is  defined  in  terms  of  the  z-component  of  the  total  wavenumber 
k.  For  a  plane  wave  having  wavenumber  components  kx  and 
kz.  Yz  is  the  ratio  of  the  z-component  of  the  particle 
velocity  to  the  pressure  (i.e..  the  longitudinal  admittance). 
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Frequency  Response  of  the  Reflection 
Coefficient  as  a  Standard  Fourier  Transform 

z 

•  z*  *  (2/kojy  kj(s)ds 

0 

z*(H) 

•  R(f,0)  *  -/  g<z*)exp[27rlfz*/c0]dz* 

0 

Time  Response  of  the  Reflection  Coefficient 
determines  the  Relative  Variation  of  Admittance 

oo 

•  •'(t.0)  -f  R(f,#)exp(-2Tift)df 

— oo 

.  _®o_  *YZ 

“  “  2YZ  d?* 

z*  =  c0t 


Viewgraph  4 

By  a  change  of  variable,  the  reflection  coefficient  can  be 
written  as  a  standard  Fourier  transform  in  the  new  depth 
coordinate  z*.  If  we  take  the  Fourier  transform  with  respect 
to  time  in  order  to  obtain  the  impulse  response  of  the  re¬ 
flection  coefficient,  we  arrive  at  the  desired  result  that  the 
impulse  response  at  time,  t,  is  related  to  the  relative  vari¬ 
ation  of  longitudinal  admittance,  Yz,  at  depth  z*  •  c0t. 
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Algorithm  to  Recover  Admittance  Profile  from  a  Single 

Impulse  Response 

•  dY2/dz*  =  -<2/Co)r  Yz,  Y^O)  =  sin  Oj  (p0c0) 

•  dz/dz*  =  (2  pc0YJ-\  z(0)  =  0 


Viewgraph  5 

These  equations  represent  Candel's,  et  al..  non-iterative 
algorithm  for  determining  admittance  as  a  function  of  depth 
from  a  single  impulse  response.  All  that  is  required  is  the 
integration  of  a  pair  of  first-order  differential  equations 
with  appropriate  initial  conditions.  Either  p(z)  or  c(z)  is 
required  to  be  known  in  order  to  effect  the  inversion.  Candel. 
et  al..  have  shown  that  if  two  impulse  responses  corresponding 
to  different  grazing  angles  are  available,  then  unique  deter¬ 
mination  of  both  p(z)  and  c(z)  is  possible  without  any  prior 
knowledge  of  either  bottom  profile.  In  this  case,  four  dif¬ 
ferential  equations  have  to  integrated  instead  of  two. 
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The  simulation  of  noisy  impulse  responses  for  the  inversion 
algorithm  is  shown  in  this  block  diagram.  For  a  given 

geoacoustic  model  and  incident  grazing  angle,  the  complex  fre¬ 
quency  response  R(f)  is  computed  at  256  discrete  equispaced 
frequencies  ranging  from  0  to  127.5  Hz.  The  sequence  is 
extended  to  1024  points  by  appending  zero  values.  A  fast 
Fourier  transform  (FFT)  is  used  to  produce  a  bandlimited 
impulse  response  r(t).  This  constitutes  a  solution  to  the 
forward  problem. 

Bach  impulse  response  is  convolved  with  a  low-pass  digital 
filter  having  a  cut-off  frequency  of  63.5  Hz.  This  is  found 

necessary  to  inhibit  Gibbs  oscillations  due  to  the  truncation 
of  the  sampled  frequency  response.  Zero-mean,  Gaussian  noise 
is  added  to  each  point  of  the  filtered  time  sequence.  Two 
"noisy"  impulse  responses,  representing  grazing  angles  60° 
and  90°,  are  supplied  to  the  inversion  algorithm. 
Reconstructed  estimates  of  the  normalized  density  p(z)//>0 

and  sound-speed  cQ/c(z)  profiles  are  compared  to  the  known 

values  comprising  the  geoacoustic  model. 
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Hatteras  Abyssal  Plain  Geoacoustic  Model 
(28°  30 '  N,  70°  30 '  W) 


Viewgraph  7 

The  geoacoustic  model  used  in  the  simulations  is  based  on 
the  interpretation  of  Herstein.  Dullea  and  Santaniello  (1979) 
for  measurements  made  at  a  site  in  the  Hatteras  Abyssal  Plain. 
Both  density  and  sound-speed  profiles  exhibit  discontinuities 
at  layer  boundaries.  The  upper  layer  is  thin  and  iso-velocity. 
The  deep  layer  is  characterized  by  an  upward  refracting 
sound-speed  gradient. 
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Definition  of  Signal-to-Noise  Ratio  (SNR) 

•  SNR  *  10  log10  (gg  dB 

ti  +  T 

•  I.  «  (1AT)  /  r2  (t)dt 

*1 

ti+T 

•  >n  =  (1/-r)  /  n2  (»)dt  %  an2 

•  SNR  *  10  iog-ioOg/cr,,2) 

GIVEN  SNR,  THEN  on  =  ls**10~SNR/2° 


Viewgraph  8 

,  this  simulation.  the  signal-to-noise  ratio  VSNR)  in  dB 

is  defined  as  10  iog10  of  the  signal  intensity  divided  by  the 
noise  intensity.  The  signal  intensity  is  computed  by  inte- 

lonroS?!^6  °-f  „the  faction  impulse  response  over  an 

appropriate  time  window  and  then  dividing  by  the  window 

length.  The  noise  intensity  is  computed  similarly  over  the 
noise  series.  Note  that  the  intensity  of  the  noise  can  be 

lt<?  variance-  Hence-  if  the  signal  intensity 
is  known,  the  variance  of  the  noise  required  to  simulate  a 
desired  SNR  can  be  computed  as  shown  in  the  bottom  line. 


Simulated  Impulse  Responses 
Grazing  Angle  0  =  60° 
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Viewgraph  9 


This  viewgraph  shows  examples  of  bandlimited  ceflectivity 
impulse  responses  generated  for  the  Hatteras  model  for  a  graz¬ 
ing  angle  of  60°.  The  upper  trace  corresponds  to  the  noise- 
free  case.  The  lower  trace  corresponds,  to  a  signal-to-noise 
ratio  of  20  dB  determined  for  an  appropriate  time  window.  T. 
The  three  "pulses"  observed  in  each  response  are  associated 
with  discontinuities  in  the  density  and/or  sound-speed  profiles 
at  layer  boundaries.  Between  the  second  and  third  "pulses" 
there  appear  continuous  reflections  arising  from  the  positive 
sound-speed  gradient  of  the  deep  layer.  Similar  responses  are 
obtained  for  a  grazing  angle  of  90°. 


9 


TD  6947 

Reconstructed  Sound-speed  Profile  Using 
Grazing  Angles  60°  and  90° 


DEPTH  <m> 


Viewgraph  10 

Using  bandlimited  impulse  responses  for  grazing  angles 
60°  and  90°.  Candel's.  et  al.,  inversion  algorithm  was  used 
to  obtain  the  reconstructed  sound-speed  profiles  shown  here. 
The  sound-speed  versus  depth  variations  are  presented  in  the 
form  of  acoustic  refractive  index  profiles.  For  comparison,  the 
Hatteras  Abyssal  Plain  input  model  values  are  also  indicated. 
The  top  comparison  relates  to  the  noise-free  case.  While  the 
global  agreement  is  excellent,  the  highest  frequency  present  in 
the  low-pass  impulse  responses  is  insufficient  to  resolve  the 
refractive  index  within  the  thin  upper  layer  of  the  model.  The 
lower  comparison  shows  the  effect  of  adding  noise  at  the  20  dB 
SNR  level.  Again,  the  global  agreement  is  good  but  the  local 
effects  of  the  noise  are  evident.  Note  that  the  sound-speed 
gradient  of  the  deep  layer  is  accurately  recovered  in  both 
cases . 
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Reconstructed  Density  Profile  Using 
Grazing  Angles  60°  and  90° 


DEPTH  (m) 


Viewgraph  11 

The  reconstructed  density  profiles  normalized  to  the 
density  in  the  region  z  <  0  are  displayed  next.  As  before,  the 
top  comparison  is  for  the  noise-free  case  and  the  lower  com¬ 
parison  is  for  a  20  dB  signal-to-noise  ratio.  These  results 
show  that  the  resolution  of  the  density  profile  is  better  than 
that  obtained  for  the  refractive  index  profile.  This  improve¬ 
ment  in  resolution  is  due  in  part  to  the  larger  relative 
variation  in  density  for  the  Hatteras  model.  Observe  that  even 
the  density  within  the  thin  upper  layer  of  the  model  is 
resolved. 
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Root  Mean  Square  Error 


Viewgraph  12 

A  set  of  simulations  was  carried  out  as  a  function  of 
signal-to-noise  ratio.  The  effect  of  signal-to-noise  ratio  on 
the  profiles  reconstructed  via  Candel's.  et  al.#  inversion 
algorithm  are  presented  in  terns  of  a  root  mean  sguare  error 
measure.  Bach  error  value  was  determined  by  summing  the 
squares  of  the  relative  differences  between  reconstructed  and 
true  profiles  over  the  depth  interval  displayed  in  the  previous 
Viewgraphs.  Error  measures  are  plotted  for  both  normalized 
density  and  refractive  index. 

For  signal-to-noise  ratios  greater  than  approximately  20  dB 
the  errors  differ  only  slightly  from  the  error  associated  with 
the  noise-free  case.  This  nonzero  error  for  the  noise-free 
case  arises  since,  first,  the  reconstructions  are  based  on 
bandlimited  impulse  responses  and,  second,  the  inversion  algo¬ 
rithm  itself  is  based  on  a  forward  scattering  approximation. 
For  signal-to-noise  ratios  less  than  IS  dB,  the  errors  increase 
rapidly.  Note  that  the  inversion  remained  stable  (i.e., 
realizable)  even  for  a  signal-to-noise  ratio  of  0  dB.  These 
results  agree  well  with  the  results  published  recently  by 
Schwetlick  (1983)  using  other  methods  of  profile  inversion. 
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Signal-to-noise  ratios  of  is  dB  have  routinely  been 
obtained  for  measurements  of  acoustic  reflectivity  at  sea.  This 
is  particularly  true  for  the  high  grazing  angles,  i.e..  close 
ranges,  appropriate  for  the  implementation  of  Candel's,  et  al.. 
inversion  method.  Replicate  averaging  or  beamforming  methods 
could  be  used  to  obtain  further  increases  in  signal-to-noise 
ratio. 

We  conclude  that  Candel's.  et  al..  non-iterative  algorithm 
is  a  viable  method  for  geoacoustic  models  similar  to  the 
Hatteras  Abyssal  Plain  model,  provided  adequate  signal-to-noise 
ratios  are  available.  An  important  feature  of  the  algorithm  is 
the  recovery  of  both  density  and  sound-speed  profiles  from  just 
two  bandlimited  impulse  responses.  No  prior  knowledge  of 
either  profile  is  required  for  the  inversion.  Beyond  this 
simulation,  we  intend  to  apply  Candel's.  et  al..  method  to  high 
signal-to-noise  measurements  of  sea-bed  reflection  responses. 


Thank  you.  Are  there  any  questions? 
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NAVSEASYSCOM  (SEA-63R;  L.  Kneesi.-63D: 

D.  Porter. -63R;  F. Romano. -63R;  R.  Farwell.-63R)  5 

NAVPGSCOL  1 

DHTNSRDC  1 

NORDA  (R.  Lauer,  320:  S.  Marshall.  115:  R.  Martin.  110A: 

L.  Solomon,  500;  R.  Gardiner.  520;  E.  Chiaka, 

B.  Bluaenthal ,  H.  Horsley.  R.  Hheatly.  530:  J. Matthews: 

G.  Stanford;  Library)  12 

NOSC  (R.  Bolam.  7133;  R.  Albrecht.  7134;  M.  Pederson. 

D.  Gordon.  712;  J.R.  McCarthy.  713.  C.  Persons.  7133; 

S.  Sullivan,  1604;  R.  Smith;  J.  Lovett.  G.  Tunstall; 

Library)  11 

NADC  (J.  Howard;  B.  Steinberg.  P.  Hass;  Library)  4 

NCSC  1 

NSWC  (R.  Stevenson;  M.  Stripling:  M.  Stallard;  Library)  4 

NRL  (O.  Diachok.  R.  Dicus,  5160;  J.  Munson.  5100. 

H.  Mosley.  5120;  A.  Eller.  5109;  Library)  6 

MPL  (V.C.  Anderson;  F.  Fischer:  B.  Williams)  3 

NAVAIR  (E.  David.  370B;  H.  Parrigian.  370J)  2 

NISC  (H.  Foxwell)  1 

DTIC  2 

ONR  (R.  Winokur.  ONR-102B;  CAPT  E.  Craig.  T.  Warfield. 

0NR-220 :  G.  Hamilton,  ONR-420;  J.  McKisic.  ONR-4250A; 

P.  Rogers.  ONR-425UA;  J.  A.  Smith:  onr-ioo;  onr-480)  9 

DIA  1 

CHESNAVFACENGCOM  1 

ARL  (Univ.  of  Texas)  (P.  Vidmar;  K.  Hawker; 

R.  Koch;  Library  4 

ARL  (Penn  State  Univ.)  (S.  McDaniel;  D.  McCammon; 

Library)  3 

Courant  Institute  (D.C  Stickler;  Library)  2 

SAI  (C.W.  Spofford.  R.  Greene)  2 

Cornell  University  (H.  Schwetlick)  1 

Northwestern  University  (G.  Kriegsmann)  1 

Univ.  of  Denver  (J.A.  DeSanto.  F.  Hagin)  2 

General  Instrument  Co.  (R.  Breton)  1 

FWG  (D.  Wills)  1 


INITIAL  DISTRIBUTION  LIST  (Coat'd) 

No.  Of 


Addressee  Copies 

CANADA 

DREP  20 

DSIS  (Microfiche  Section,  Report  Collection)  3 

DREA  1 

ORAE  Library  (DMOR)  2 
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CDLS/W  CDR  1 

CMDO  1 

DMRS  1 

Commander  Maritime  Command  (MC/ORD,  SSO  Ocean)  2 

CFMWS  1 

Maritime  Headquarters  Pacific  (SSO  Op  Rsch)  1 

DST(SE) 3  1 

DTA(M) 3  1 

RRMC  (Dept,  of  Oceanography)  1 

D  Met  Oc  1 

Project  Officer  IEP  ABCANZ-2.  NOP  981.  Pentagon  3 

BRITAIN 

DRIC  3 

AUWE  1 

RAE  1 

AUSTRALIA 

DRC  1 

NEW  ZEALAND 

DSE  1 


